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r^ ■ Abstract 

O 

o: 

"^ , Lifetime and oscillations of B^ mesons have been studied in events with a large 

D I transverse momentum lepton and a Dg of opposite electric charge in the same 

p ■ hemisphere, selected from about 3.6 million hadronic Z° decays accumulated 

O : by DELPHI between 1992 and 1995. 

The Bg lifetime and the fractional width difference between the two physical 

B° states have been found to be: 



H ; tbo = {lA2toltistat.) ± Om{syst.)) ps 

APbo/Pbo < 0.46 at the 95% C.L. 

In the latter result it has been assumed that tro = Tro . 

Using the same sample, a limit on the mass difference between the physical Bg 
states has been set: 

AmRO > 7.4 ps"^ at the 95% C.L. 

with a corresponding sensitivity equal to 8.1 ps^^. 



(Eur. Phys. J. C16(2000)555) 
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1 Introduction 

In this paper, the average hfetime of the Bg meson has been measured and hmits have 
been derived on the oscillation frequency of the Bg-B^ system, Am^o, and on the decay 
width difference, AFbo, between mass eigenstates of this system. 

Starting with a B° meson produced at time t=0, the probability, V, to observe a Bg or a 
B^ decaying at the proper time t can be written, neglecting effects from CP violation: 

^[B° - B°(BO)] = ^e-^-2'[cosH^t) ± cos(AmBot)] (1) 

where Fbo = (Fro + Tko)/2, AFbo = F^o — F^o and ArriBO = mUo — mko. L and 
H denote the light and heavy physical states, respectively; AFbo and AmBO are defined 
to be positive [^ and the plus (minus) signs refer to B^ (B°) decays. The oscillation 
period gives a direct measurement of the mass difference between the two physical states. 
The Standard Model predicts that AFbo ^ Am^o, for which the previous expression 
simplifies to : 

p™.. _ p(bO^BO) = FBoe-^Bot^os2(^^^) (2) 

and similarly: 

V^- =P(B°^BO) = FBoe-^Bot,i^2(AmBet^ ^^^ 

The oscillation frequency, proportional to Am^o, can be obtained from the fit of the time 
distributions given in relations (||) and (0), whereas expression ([^), without distinguish- 
ing between the B^ and the B°, can be used to determine the average lifetime and the 
difference between the lifetimes of the heavy and light mass eigenstates. 

B physics allows a precise determination of some of the parameters of the Cabibbo 
Kobayashi Maskawa (CKM) matrix. All the nine elements can be expressed in term 
of four parameters that are, in Wolfenstein parametrization 0, A, A, p and t]. The val- 
ues of p and 7] are the most uncertain. 

Several quantities which depend on p and r] can be measured and, if the Standard Model 
is correct, they must define compatible values for the two parameters, inside measurement 
errors and theoretical uncertainties. 

These quantities are ex, the parameter introduced to measure CP violation in the K 
system, |Vufe|/|Vcfe|, the ratio between the modulus of the CKM matrix elements corre- 
sponding to 6 — s> M and b ^ c transitions and the mass difference Atji-qo . 

In the Standard Model, B^-B^ {q = d, s) mixing is a direct consequence of second order 
weak interactions. Having kept only the dominant top quark contribution, Attt-bo can be 
expressed in terms of Standard Model parameters 0: 

^^Bo = ^,\Vt>,m,\'m^,m^JlB^^VBFi^). (4) 

2 

In this expression Gp is the Fermi coupling constant; F(xt), with Xt = -^, results from 
the evaluation of the second order weak "box" diagram responsible for the mixing and 
has a smooth dependence on Xt] tjb is a QCD correction factor obtained at next to leading 
order in perturbative QCD |Q. The dominant uncertainties in Equation (^) come from 
the evaluation of the B meson decay constant /b and of the "bag" parameter Bb ■ 



The mass differences Am^o and Atjibo involve the CKM elements Vtd and Vts- Neglecting 

d ^ 

terms of order A^, these are given by: 



\Vtd\ = AX'^{l-py + r]^ ; \Vts\ = Ay. (5) 

In the Wolfenstein parametrization, \Vts\ is independent of p and rj. A measurement 
of Am-BO is thus a way to measure the value of the non perturbative QCD parameters. 
Direct information on Vtd can be inferred by measuring Am^o . 

Several experiments have accurately measured Am-^o , nevertheless this precision cannot 
be fully exploited to extract information on p and rj because of the large uncertainty 
which originates in the evaluation of the non-perturbative QCD parameters. 
An efficient constraint is the ratio between the Standard Model expectations for Am^o 
and Att^bo, given by: 

AmBO "^Bo/lo^Bor/BO l-^ |2 

AmBO mBo/lo^Bor^BO \Vtsf 

A measurement of the ratio ATn^o/AmBO gives the same type of constraint, in the p — r] 
plane, as a measurement of Am^o, but because only ratio /bo//bo and B^o/B^o are in- 

d d s d s 

volved, some of the theoretical uncertainties cancel . 

Using existing measurements which constrain p and rj, except those on AmBO, the dis- 
tribution for the expected values of Attibo can be obtained. It has been shown, in the 
context of Standard Model and QCD assumptions, that Att^bo has to lie, at 68% C.L., 
between 12 and 17.6 ps~^ and is expected to be smaller than 20ps~^ at 95% C.L. @|. 

The Bg meson lifetime is expected to be equal to the B^ lifetime [|l| within one percent. 
In the Standard Model, the ratio between the mass difference and decay width in the 
6*^-6° system is of the order [mb/mty, although large QCD corrections are expected. 
Explicit calculations to leading order in QCD correction, in the HQE (Heavy Quark 
Expansion) formalism |T|, predict: 

ArBo/rBO = o.i6+°:S^ 

where the quoted error is dominated by the uncertainty related to hadronic matrix ele- 
ments. 
Recent calculations [0 at next-to-leading order predict a lower value: 



AEBo/rBO = 0.054+0; 



0.016 
032 



An interesting approach consists in using the ratio between AFbo and Am^o 0: 

(2.63l?iDlO-' (7) 



^^60 ,^ ^^+o.67^ 



Attibo 

to constrain the upper part of the Am-BO spectrum with an upper limit on AFbo/Fbo. 
If, in future, the theoretical uncertainty can be reduced, this method can give an alterna- 
tive approach in determining Attt-bo via AFbo and, in conjunction with the determination 
of Am^o , can provide an extra constraint on the p and r] parameters. 

The results presented in the following have been obtained from data accumulated by 
DELPHI experiment at LEP between 1992 and 1995, corresponding to about 3.6 million 
hadronic Z^ decays. The main features of these analyses are: 



• a precise measurement of the B decay proper time; 

• a determination of the charge of the b quark at the B- meson decay time (decay tag); 

• a determination of the sign of the b quark at production time (production tag). 

The first item is common to the three studies on AmBO, tbo and AFbo while the others 
are specific to the oscillation analyses. For these last, the principle of the measurement 
is as follows. Each of the charged and neutral particles measured in an event is assigned 
to one of the two hemispheres defined by the plane transverse to the sphericity axis. A 
"production tag" is used to estimate the b/b sign of the initial quark at the production 
point. The decay time of the B hadron is evaluated and a "decay tag" is defined, correlated 
with the b/b content of the decaying hadron. The analysis is performed using events 
containing a lepton emitted at large transverse momentum, pt, relative to its jet axis 
accompanied by an exclusively (or partially) reconstructed Dg in the same hemisphere 
and of opposite electric charge. The lepton charge defines the "decay tag". Different 
variables defined in the same and in the opposite hemisphere, are used to determine the 
"production tag" . 
Similar analyses have been performed by the ALEPH, CDF and OPAL Collaborations 

mm- 

Section ^ describes the main features of the DELPHI detector, the event selection and the 
event simulation. Section |^ describes the selection of the Dgi sample. Section ^ presents 
the Bg lifetime measurement. Section ^ presents the result on the lifetime difference. 
Section ^ is devoted to the study of Bg-B° oscillations with the Dgi sample: the first 
part describes the "production tag" algorithm while the second part presents the fitting 
procedure and the result on Atjibo. 

2 The DELPHI detector 

The events used in this analysis have been recorded with the DELPHI detector at LEP 
operating at energies close to the Z^ peak. The DELPHI detector and its performance 
have been described in detail elsewhere fill]. In this section are summarized the most 
relevant characteristics for this analysis. 

2.1 Global event reconstruction 

2.1.1 Charged particles reconstruction 

The detector elements used for tracking are the Vertex Detector (VD), the Inner 
Detector (ID), the Time Projection Chamber (TPC) and the Outer Detector (OD). 

The VD provided the high precision needed near the primary vertex. For the data 
taken from 1991 to 1993, the VD consisted of three cylindrical layers of silicon detectors 
(radii 6.3, 9.0 and 10.9 cm) measuring points in the plane transverse to the beam direction 
{r(f) coordinate) in the polar angle range 43° < 9 < 137°. In 1994, two layers have 
been equipped with detector modules with double sided readout, providing a single hit 
precision of 7.6 /zm in the r0 coordinate, similar to that obtained previously, and 9 /im in 
the coordinate parallel to the beam {z) []T2| . For high momentum particles with associated 
hits in the VD, the extrapolation precision close to the interaction region is 20 /xm in the 
r0 plane and 34 /xm in the rz plane. 

Charged particle tracks have been reconstructed with 95% efficiency and with a 
momentum resolution Op/p < 2.0 x 10~^p [p in GeV/c) in the polar angle region 
25° < ^ < 155°. 



2.1.2 Energy reconstruction 

The total energy in the event is determined by using all information available from the 
tracking detectors and the calorimeters. For charged particles, the momentum measured 
in the tracking detector is used. Photons are detected and their energy measured in the 
electromagnetic calorimeters, whereas the hadron calorimeter detects long lived neutral 
hadrons such as neutrons and -ft'^'s. 

The electromagnetic calorimetry system of DELPHI is composed of a barrel calorime- 
ter, the HPC, covering the polar angle region 46° < 6 < 134°, and a forward calorimeter, 
the FEMC, for polar angles 8° < 9 < 35° and 145° < 9 < 172°. The relative precision on 
the measured energy E has been parametrized as cte/E = 0.32/ y/E (B 0.043 {E in GeV) 
in the barrel, and (Te/ E = 0.12/vE (B 0.03 {E in GeV) in the forward region. 

The hadronic calorimeter, HCAL, has been installed in the return yoke of the DELPHI 
solenoid. In the barrel region, the energy has been reconstructed with a precision of 
aE/E = 1.12/ VE © 0.21 {E in GeV). 

2.1.3 Hadronic Z° selection 

Hadronic events from Z^ decays have been selected by requiring a charged multiplicity 
greater than four and a total energy of charged particles greater than 0.12-^5, where -y/i 
is the centre-of-mass energy and all particles have been assumed to be pions; charged par- 
ticles have been required to have a momentum greater than 0.4 GeV/c and a polar angle 
between 20° and 160°. The overall trigger and selection efficiency is (95.0±0.1)% |jl3|- A 
total of about 3.6 million hadronic events has been obtained from the 1992-1995 data. 

2.2 Particle identification 

2.2.1 Lepton identification 

Lepton identification in the DELPHI detector is based on the barrel electromagnetic 
calorimeter and the muon chambers. Only particles with momentum larger than 2 GeV/c 
have been considered as possible lepton candidates. 

Muon chambers consisted, in the barrel region, of three layers covering the polar 
regions 53° < 9 < 88.5° and 91.5° < 9 < 127°. The first layer contained three planes 
of chambers and was inside the return yoke of the magnet after 90 cm of iron, while 
the other two, with two chamber planes each, were mounted outside the yoke behind a 
further 20 cm of iron. In the end-caps there were two layers of muon chambers mounted 
one outside and one inside the return yoke of the magnet. Each consisted of two planes 
of active chambers covering the polar angle regions 20° < 9 < 42° and 138° < 9 < 160° 
where the charged particle tracking was efficient. 

The probability of a particle being a muon has been calculated from a global x^ of 
the match between the track extrapolation to the muon chambers and the hits observed 
there. Four identification flags are given as output of the muon identiflcation in decreasing 
order of efficiency: very loose, loose, standard and tight. In this analysis the loose 
selection has been applied corresponding to an efficiency of (94.8 ±0.1)% with a hadron 
misidentification probability of (1.5 ± 0.1)%. 

Electron identification has been performed using two independent and complementary 
measurements, the dE/dx measurement of the TPC (described in Section 2.2.2) and the 
energy deposition in the HPC. Probabilities from calorimetric measurements and tracking 
are combined to produce an overall probability for the electron hypothesis. Three levels 



of identification are given: loose, standard and tight. 

The loose selection has been applied for this analysis corresponding to an efficiency of 

80 % with an hadron misidentifation probability of ~ 1.6 %. 

2.2.2 Hadron identification 

Hadron identification relied on the RICH detector and on the specific ionization mea- 
surement performed by the TPC 

The RICH detector |jT^ used two radiators. A gas radiator separated kaons from 
pions between 3 and 9 GeV/c, where kaons gave no Cherenkov light whereas pions did, 
and between 9 and 16 GeV/c, using the measured Cherenkov angle. It also provided 
kaon/proton separation from 8 to 20 GeV/c. A liquid radiator, which has been fully 
operational for 1994 and 1995 data, provided p/K/t^ separation in the momentum range 
1.5-7 GeV/c. 

The specific energy loss per unit length (dE/dx) is measured in the TPC by using up 
to 192 sense wires. At least 30 contributing measurements have been required to compute 
the truncated mean. In the momentum range 3 < p < 25 GeV/c, this is fulfilled for 55% 
of the tracks, and the dE/dx measurement has a precision of ±7%. 

The combination of the two measurements, dE/dx and RICH angles, provides three 
levels of pion, kaon and proton tag (loose, standard, tight) corresponding to different 
purities. A tag for "Heavy Particle" is also given in order to separate pions from heavier 
hadrons with high efficiency. 

The Standard "Heavy Particle" fiag has an efficiency of about 70 % with a pion misidenti- 
fication probability of 10 % for charged particle with momentum greater than 0.7 GeV/c. 

2.2.3 A° and K^ reconstruction 

The A'^ -^ pTT~ and K^ — > tt+tt^ decays have been reconstructed if the distance in the 
r0 plane between the V^ decay point and the primary vertex is less than 90 cm. This 
condition meant that the decay products have track segments at least 20 cm long in the 
TPC. The reconstruction of the V° vertex and selection cuts are described in detail in 



reference 1 1 



Only K candidates passing the "tight" selection criteria have been retained for this 
analysis. 

2.2.4 7r° reconstruction 

The 77° — > 77 decays are reconstructed by fitting all 77 pairs whose invariant mass is 
within 20 MeV of the nominal vr" mass, using the nominal n^ mass as a constraint. The 
fit probability has to be larger than 1%. 

2.3 Primary vertex reconstruction and event topology 

The location of the e^e~ interaction has been reconstructed on an event-by-event basis 



using the beam spot position as a constraint [|IT]. In 1994 and 1995 data, the position 
of the primary vertex transverse to the beam has been determined with a precision of 
about 40 /^m in the horizontal direction, and about 10 /im in the vertical direction. For 
1992 and 1993 data, the uncertainties are larger by about 50%. 

Each selected event has been divided into two hemispheres separated by the plane 
transverse to the sphericity axis. A clustering analysis based on the JETSET algorithm 



LUCLUS |T5[ with default parameters has been used to define the jets, using both charged 



and neutral particles. These jets have been used to measure the Pt of each particle in the 
event, defined as its momentum transverse to the axis of the rest of the jet it belonged 
to, after removing the particle itself. 

The different detector configurations, both for hadron identification and vertex reso- 
lution, implies, in the rest of the analysis, a separate treatment of the data taken before 
and after 1994. 

2.4 6-tagging 

The 6-tagging package developed by the DELPHI collaboration has been described in 



reference ||T6[. The impact parameters of the charged particle tracks, with respect to the 
primary vertex, have been used to build the probability that all tracks come from this 
vertex. Due to the long B-hadrons lifetime, the probability distribution is peaked at zero 
for events which contained beauty whereas it is fiat for events containing light quarks. 
The 6-tagging algorithm has been used in this analysis to select control samples with low 
b purity. 

2.5 Event simulation 



Simulated events have been generated using the JETSET 7.3 program |T5[ with pa- 



rameters tuned as in Ijl^ and using an updated description of B decays. B hadron 



semileptonic decays have been simulated using the ISGW model [1^]. Generated events 



have been followed through the full simulation of the DELPHI detector (DELSIM) |Tl| 
and the resulting simulated raw data have been processed through the same reconstruc- 
tion and analysis programs as the real data. 



3 The D=^£^ sample selection 



B° semileptonic decay^ have been selected requiring the presence of a D+ meson 
correlated with a high pj- lepton of opposite electric charge in the same hemisphere: 

BO -. Dtrvix. 

The Ds mesons have been reconstructed in six non-leptonic and two semileptonic decay 
channels: 

D+^(t)n+ 0^K+K-; 

D3+ -^ T^K+ ir° -^ K-7r+; 

D+ -^ KgK+ Kg -^ TT+TT-; 

D3+ ^ </)7r+7r-7r+ ^ K+K"; 

D+^(t)n+TT° 0^K+K-; 

D+^K''°K^+ ir°^K-7r+, K*+^Kg7r+; 



D3+ -^ (Pe+Ue (p -^ K+K- 
D+ -^ 0/i+z/^ ^ K+K- 



^Charge conjugation is always implied. 



In addition, partially reconstructed D^ have been selected requiring the presence of a 
meson (reconstructed in the K+K" decay channel) accompanied by an hadron /i"*" in the 
same hemisphere. 



DJ 



(ph+X 



In the following the first eight decay modes will be referred as the Ds 
last one as the (pih sample. 



sample and the 



3.1 Selection of the 



J7T 



K*\- 



KgK- 



and (j)£^iy decay modes 



Each Ds decay mode has been reconstructed by making all possible combinations of 
particles in the same hemisphere. In D+ semileptonic decays, the ambiguity between the 
two leptons has been removed by assigning the lepton to the D^^ (B^) if the mass of the 
(pi system, M{(j)i), is below (above) the nominal D^ mass. If the two leptons both gave 
a M{(j)i) above or below the Dg mass, the event was rejected. 

The measured position of the D+ decay vertex and momentum together with their mea- 
surement errors, have been used to form a new track (called pseudo-track) that contains 
the measured parameters of the D+ particle. 

A candidate B° decay vertex has been obtained by intercepting the D+ pseudo-track with 
the one of a lepton. To guarantee a precise determination of the position of this secondary 
vertex, at least one VD hit has been required to be associated to the lepton and to at 
least two tracks from the D+ decay products. The x^ of the reconstructed D^ and B^ 
vertices have been required to be smaller than 40 and 20 respectively. 
In order to suppress fake leptons and B hadron cascade decays (6 — >■ c — >■ i~^), additional 
selection criteria have been applied to the Dgi pairs, which are summarized in Table |^. 
For the channel D+ — > (pi'^iy requirements on the (pii mass and momentum have been re- 
duced as compared to the other channels to account for the additional escaping neutrino. 
Due to the smaller combinatorial background under the Dg signal, in the Dg -^ (pn^ and 
Dg — > 0^+z/ decay channels, the pt cut has been lowered to 1 GeV/c. 





07r+ 


0£+ 


Others 


PT{e){GeV/c) 

M{BJ){GeV/c^) 

P{DJ){GeV/c) 


> 1 

G [3,5.5] 
>14 


> 1 

G [2.5,5.5] 
> 12 


> 1.2 

G [3,5.5] 

> 14 



Table 1: Selection criteria applied to the lepton and Dg candidates. 



A tighter selection was then applied, separately for each decay mode, using a discrim- 
inant variable built with the variables listed in Table 0. 
These variables are: 

• the momenta, P, and masses, M, of the decay products; 

• the cosine of the helicity angle, ip, for the (l)7i~^ and K K+ decay modes; 

• Hj£,, defining whether the hadron identification from Section 2.2.2| favours the vr, K 
or proton hypothesis; 

• LjD, defining whether the lepton identification from Section |2.2.1| identifies a particle 
from the D+ semileptonic decay as an electron or a muon (used only for leptons 
coming from the D+ semileptonic decays). 



For each quantity the probabihty densities for the signal (S) (Dg^ from B^ seniileptonic 
decays) and for the combinatorial background (B) (fake D^i candidates in qq events) have 
been parametrized using the simulation; the discriminant variable Xd^ is then defined as 



^=n^ =n 



bi[Xi) 

Bi{xi) 



Xr 



R 



R+l 



where i runs over the number of variables (which actual values are Xi). The combinatorial 
background is concentrated close to Xd, = while the Dg signal accumulates close to 
Xd, = 1. The definition of Xd, provides an optimal separation between the signal and 
the combinatorial background if the individual discriminant variables Xi are independent; 
in case of correlations the separation power decreases but no bias is introduced. 



(j)TX 


KO*K 


K°sK 


(j)£+ 


pm 


^(Ds) 


^(Ds) 


P{4>) 


pi<p)/pm 


P{K*')/P{D,) 


P(KO)/P(D,) 




HidK, 


HidK^ 


HidK 


HidK, 


Hid K2 


Hid K2 




Hid K2 


Hid tt 


Hid vt 




Lid im 


cos (■?/') 


cos{%l>) 






M(0) 


M(K*0) 







Table 2: List of the quantities which are used, in the different decay channels, to construct 
a discriminant variable between Bg semileptonic decays and background events. 



The distributions of this variable obtained in data and in the simulation are shown in 
Figure |I| for the ^vr"*" decay channel. 

The optimal value of the cut on the discriminant variable has been studied on simulated 
events, separately for each channel and for each detector configuration, in order to keep 
high efficiency (Table ^. A very loose cut has been applied on the ^vr"*" channel because 
of its small combinatorial background. 

The individual event purity has been evaluated, in the following, from the distribution 
of the discriminant variable for signal and combinatorial background. 





07r+ 


KO*K 


KgK 


0£+ 


92-93 
94-95 


>0.05 
>0.03 


>0.75 
>0.85 


>0.80 
>0.90 


>0.75 
>0.90 



Table 3: Values of the cuts applied on the discriminant variable X^^ to select Bg semilep- 
tonic decay candidates. 



In addition, for the two channels (K*''K and KgK), which receive contributions from 
kinematic reflections of non strange B decays, the bachelor kaon has been required to be 
incompatible with the pion hypothesis. 
Further background suppression has been obtained by placing a requirement on the Dg 



flight distance L(Ds). The small effect induced on the decay time acceptance has been 
taken into account in the following. This requirement has been applied, depending on the 
resolution on the decay distance observed in the different Dg decay channels and on the 
level of the combinatorial background: L(I)s) > for 07r and K**'K"'", L(Ds)/o"(L(Ds)) > 
-3 for KgK+ and LiD,)/a{L{D,)) > -1 for (j)i+. 

Finally, for the semileptonic decay modes (with two neutrinos in the final state) an 
algorithm has been developed to estimate the missing energy, Emiss, defined as: 

where the visible energy {Ems) is the sum of the energies of charged particles and photons 
in the same hemisphere as the Dgf candidate. Using four-momentum conservation, the 
total energy {Etot) in that hemisphere is: 

■'tot -^beam 



T-i j-i I same opp 

^tot — -^beam + 



^-^beam 

where Msame and Mopp are the hemisphere invariant masses of the same and opposite 
hemispheres respectively. A positive missing energy Emiss has been required. 

3.2 Selection of the ^tt+tt+tt , ^tt+tt^ and K K*+ decay modes 

These three decay modes have been searched for in the 94 and 95 data only. 
Ds£ pairs have been selected by requiring M(Ds£) > 3.0 GeV/c^, Pt{^) > 1-2 GeV/c and 
X^ (Ds£ vertex) < 20 (except for the (j)7i^7i~^7i~ decay mode in which no x^ cut has been 
applied). 

In each event only one candidate is kept. The procedure is the following: if more than one 
candidate passed all the selection criteria only the one with the highest lepton transverse 
momentum and, if the same lepton candidate is attached to several D+ candidates the 
highest D+ momentum, is kept. 

It has been verified that this requirement keeps the signal with high efficiency and removes 
some of the combinatorial background. 

3.2.1 D+ -. K*°K*+ 

D+ candidates have been selected by reconstructing K — > K^7r+ and K*+ -^ K^vr"^ 
decays. Kg candidates have been reconstructed in the mode K^ — > tt+tt^ by combining all 
pairs of oppositely charged particles and applying the "tight" selection criteria described 
in [|1T|. The K^ has been then combined with two charged particles of the same sign, and a 
third with opposite charge. If more than one D+ candidate could be reconstructed by the 
same four particles (by swapping the two pion candidates for example) the D+ candidate 
minimizing the squared mass difference (M(K-7r+) - M(K*°))2 + (M(K°7r+) - M{K*+)f 
has been chosen, where M(K ) and M(K*"'") are the nominal K* masses [|T9|. The three 
charged particle tracks have been fitted to a common vertex and the x^ of this vertex has 
been required to be smaller than 30. To improve the resolution on the vertex position, 
all three tracks have been required to have at least one VD hit. 

K~7r"'" and Kg7r+ mass combinations have been selected if their effective masses are within 
±75 and ±95 MeV/c^ of the nominal neutral and charged K* mass respectively. 
The charged pion and kaon from K* decays must have a momentum larger than 1 and 
1.5 GeV/c respectively. The charged and neutral K* mesons must have a momentum 
larger than 4 and 3.5 GeV/c respectively and D^ mesons have a momentum larger than 
11 GeV/c. 
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3.2.2 D+ -^ (f)7T7V7T 



s 



The is reconstructed in the decay channel — >■ K+K" by taking all possible pairs 
of oppositely charged particle tracks that have an invariant mass within 13 MeV/c^ 
of the nominal meson mass |T3. Neither kaon candidate should be tagged by the 



combined RICH and dE/dX measurements as pions ("tight" selection). Three tracks, 
each compatible with the pion hypothesis as given by the combined RICH and dE/dX 
measurements, have been then added to the candidate to make a D+. The five tracks 
have been required to be compatible with a single vertex, but no requirement has been 
applied on the x^ of the vertex fit. Three of the five tracks have been required to have 
at least one VD hit and two of the three pion candidates have been required to have a 
momentum above 1.2 GeV/c. 

In addition, kaons from the decays must have a momentum larger than 1.8 GeV/c. 
Individual pion momenta must be larger than 700 MeV/c and the Dg candidate momen- 
tum must be larger than 9 GeV/c. 

3.2.3 D+ -^ 07r7r° 

The is reconstructed using the same selection criteria as for the previous channel. A 
third track, which has been required not to be tagged as a kaon by the combined RICH 
and dE/dx, and a reconstructed 7r° (Section. p.2.4|) have been added to the candidate. 



The three charged tracks have been fitted to a common vertex. To improve the resolution 

on the vertex position, each of the three tracks has been required to be associated to at 

least one VD hit each. 

In addition, kaons from the decay must have a momentum larger than 2.5 GeV/c. 

The momentum of the charged pion and of the Dg must be larger than 1 and 10 GeV/c 

respectively. 

3.3 Summary for the D^i selected events 

3.3.1 Non leptonic Dg modes 

In the D^ mass region, an excess of "right-sign" (Dfi^) over "wrong-sign" (Dfi^) 
combinations is observed in each channel (Figure 0). The estimated number of signal 
events and the yields for the combinatorial background in all the studied modes are 
summarized in Table ^. The mass distribution for non-leptonic decays has been fitted 
with two Gaussian distributions of equal widths to account for the D^ and D"*" signals 
and a polynomial function for the combinatorial background. The D"*" mass has been 
fixed at the nominal value of 1.869 GeV/c^ [|T^. The overall mass distribution for non- 
leptonic decays is shown in (Figure ^). The fit yields a signal of (206 ± 21) Dg decays in 
"right-sign" combinations, centred at a mass of (1.9680 ± 0.0016) GeV/c^ with a width 
of (14±1) MeV/c2. 

3.3.2 Semileptonic Dg modes 

Selected events show an excess of "right-sign" with respect to "wrong-sign" combi- 
nations (Figure ^). The K^K^ invariant mass distribution for "right sign" events has 
been fitted with a Breit-Wigner distribution to account for the signal and a polynomial 
function to describe the combinatorial background. The fit gives (80 ± 16) events (see 
Table |) centred at a mass of (1.020 ± 0.001) GeV/c^ with a total width (F) of (5 ± 1) 
MeV/c2. 
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Dg decay modes 


Estimated signal 


Combinatorial background / Total 


Ds -^ 07r+ 


83 ±11 


0.38 ±0.06 


D, ^ K^^K+ 


60 ±11 


0.45 ±0.06 


Ds ^ KOK+ 


22 ±7 


0.48 ±0.10 


Ds ^ ir^K*+ 


21 ±5 


0.31 ±0.07 


Ds ^ (^7r"'"7r"'"7r~ 


10 ±4 


0.39 ±0.10 


Ds -^ (pn+n^ 


18 ±6 


0.39 ±0.10 


Ds ^ 0£+z/ 


80 ±16 


0.38 ±0.06 



Table 4: Numbers of Dg signal events and fractions of combinatorial background events 
measured in the different Ds decay channels. The level of the combinatorial background 
has been evaluated inside a mass interval of±2a (^±1.5r^ centred on the measured Dg ((f)) 
mass. 



3.4 Selection of the cjylh inclusive channel 

Inclusive Bg semileptonic decays are reconstructed by requiring, in the same hemi- 
sphere, a high pt lepton and a reconstructed -^ K+K^. This analysis is expected to 
be more efficient than analyses based on completely reconstructed D+, at the cost of a 
higher background. The extra contamination comes mainly from combinatorial K+K" 
pairs and from non-strange B-decays. 

In order to avoid a statistical overlap with the Ds£ sample considered previously, all 
K+K~£ triplets selected in the Dg^ channels containing a in the final state have been 
excluded from the present sample. 

The analysis of the 0£/i channel has been performed using 94-95 data only. 
Leptons are required to have a momentum and a transverse momentum larger than 
3.0 GeV/c and 1.0 GeV/c respectively. A pair of oppositely charged identified kaons is 
considered as a candidate provided their combined momentum is above 3.0 GeV/c. 
Considering the remaining particles of charge opposite to the lepton, the hadron h with 
the highest momentum projected along the (j) direction is associated to the D+ decay ver- 
tex. The K+K~h+ vertex is fitted, and the D^ pseudo-track is reconstructed and fitted 
with the lepton track to estimate the B decay vertex. The mass distribution of the K+K" 
pairs has been fitted with a Breit-Wigner function to account for true mesons and a 
polynomial function for the combinatorial background (Figure ^. 

Accepting events within ±ir of the fitted mass, where F corresponds to the fitted 
width of the signal, 441 events are retained, including a combinatorial background of 
(45.2 ±4.5)%. 



3.5 Sample composition 

The lifetime and the oscillations of Bg mesons have been studied selecting, in the 
Dg^ sample, right-sign events lying in a mass interval of ±2cr (±1.5F) centered on the 
measured Dg (0) mass and, in the (j)lh sample, events with the candidate meson in a 
mass interval of ±1F centered on the measured mass. 
The following components, entering into the selected sample, have to be considered: 
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(p-K 


(/)£ 


Others 


fbcll fhf 


0.151 ±0.018 


0.148 ±0.025 


0.114 ±0.020 



Table 5: Ratio between DgD and signal yields in the three Dgf classes, 
fbkg'- fraction of candidates from the combinatorial background: it has been evalu- 



ated from the fit of the mass distributions on Ds^ and 



events; 



• fff fraction of candidates coming from events having a fake lepton and a real Dg 
or meson (in the (j)£h analysis this category includes also events containing true 
leptons and (j) mesons coming from charm decays or light quark hadronization) ; 

• fhcf fraction of candidates in which the high pt lepton originates from a "cascade" 
decay {b ^ c ^ £); 

• f^f fraction of semileptonic decays of non-strange B mesons 

B" n 

• /j,/ : fraction of semileptonic decays of the B^ meson. 

Only the last four components (i.e. background and signal coming from physical pro- 
cesses) will be detailed in the following: the estimation of the combinatorial background 
has been already reported in previous sections. 



3.5.1 Composition of the Dgi sample 

In the Dgi sample the Dg signal of the "right" sign correlation is dominated by B^ 
semileptonic decays; other minor sources of Dg^ candidates are: 

• //^ 

a possible contribution from this source (D+-fake £) would give the same contribution 

in right and wrong sign candidates. Since no excess has been observed in wrong sign 
candidates this component has been neglected. 

• fbcl- 

it is the expected fraction of "cascade" decays (B ^ D Dg'-^-'^X) followed by the 
semileptonic decay D -^ £~VX yielding right-sign Dfi^ pairs (referred also as /dsd)- 
This background corresponds approximately to the same number of events as the 
signal ||2^, but the selection efficiency is lower because of the requirement of a high 
Pt lepton and of a high mass of the(Ds£) system. These selection criteria reduce 
the DgD background fractions to the values reported in Table ^. Quoted errors 
on these fractions result from the uncertainties on the branching fractions of the 
contributing processes and from the errors on the respective experimental selection 
efficiencies. 

two contributions to this fraction have been considered: 

~ frefi'- the fraction of events from D+ -^ K^Tr+Tr"*" and D+ -^ Kgyr"*" decays in 
which a vr"*" has been misidentified as a K+ which give candidates in the Dg 
mass region. If the D"*" is accompanied by an oppositely charged lepton in the 
decay B^^ -^ D+£~z7X, it looks like a B^ semileptonic decay. The fractions 
frefi/fB, = 0.054 ± 0.015 and frefi/fB, = 0.069 ± 0.025 have been estimated for 
the K K"*" and KgK"*" decay channels, respectively. 

— A D^^"^ pair from a non-strange B meson decay, with the lepton emitted from 
a direct B semileptonic decay, may come from the decay B —>■ DsKX^'F. The 
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(f)7T 


(pi 


KgK+ 


K*ok+ 


Others 


/bo 


0.869 ±0.014 


0.871 ±0.019 


0.845 ± 0.023 


0.856 ±0.018 


0.898 ±0.016 


/dsD 


0.131 ±0.016 


0.129 ±0.022 


0.096 ±0.021 


0.098 ±0.020 


0.102 ±0.018 


frefl. 


- 


- 


0.059 ±0.022 


0.046 ±0.014 


- 



Table 6: Estimated composition of the Dg signal in the Dgi sample 



production of Dg in B decays not originating from W~^ —>■ cs, has been measured 
by CLEO |21], but no measurement of this production in semileptonic decays 
exists yet. This process implies the production of a D** followed by its decay 



into DgK. This decay is suppressed by phase space (the Dg K system has a large 
mass) and by the required additional ss pair. A detailed calculation shows that 
the contribution of this process is [O: 



Br(6 -^ B 



Br( 



B° 



' < 10%. 



D«£-z7) 



Assuming a selection efficiency similar to the one for the DgD component the 
contribution of this decay channel is below 2% and, for this reason, has been 
neglected in the following. 

Taking into account the above components, the estimated number of Bg semileptonic 

decays in the sample of 436 candidates is 230 ± 18. 

The signal composition for each Dg decay mode is given in Table ^. 

In order to increase the effective Bg purity of the selected sample, signal and back- 
ground fractions have been calculated on an event by event basis using the probability 
density functions oipril) and Xd^ (defined in Section p.l| ): 

feff 
J bkg 
feff 
/bo 
eff 



fejj 

Jd^d 

feff 
Jrefl 



fbkg^Comb{X-o^)Tcomb (Pt) /Tot 
/bo J^D, (Xd J J^BO {PT)/T0t 
/D,D^D,(XDjJ^DsD(pr)/Tot 
frefl^H, (XdJJ^BO {PT)/T0t 



where JF^^, Tcomb-, -^DsD, ^bo are the probability densities for the Dg mesons, the combi- 
natorial background, the DgD background and the B^ signal events, respectively. 
In these expressions. Tot is a normalisation factor such that: 



feff I feff I ceff 
Jbkg "T JbO "t J DiiD 



+ /, 



eff 
refl 



1. 



The distributions of the values of the Xd^ and pt variables are shown in Figure ^. The 
use of this procedure is equivalent to increasing the statistics by a factor 1.2. 



3.5.2 Composition of the (pih sample 



The different contributions to (pdh candidates contained in the selected K+K" mass 
interval of ±6.6 MeV/c^ around the nominal mass and corresponding to a real meson 
are shown in Table |^ and have been estimated using simulated events and measured 
branching fractions. Quoted uncertainties originate from the finite Monte Carlo statistics. 
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Source 


(%) 


ffi 


22.5 ±2.4 


J be 


49.4 ±2.1 


fbcl 


11.3 ±1.2 


Jbi 


16.9^1° 



Table 7: Estimated composition of the (p signal in the (pih sample 



except for one attached to the signal fraction which is dominated by f^o x Br(Bg 



D+rz/^X) = (0.86±0.09l°:i)% @- 



The number of B^ semileptonic decays contained in this sample has been evaluated to 



be 4llig. 



3.6 Measurement of the B meson decay time 

For each event, the B^ decay time is obtained from the measured decay length (Lbo) 



and the estimate of the B^ momentum [pbo] 
measured energies: 



The Bg momentum is estimated using the 



Pbo' = {E{DJ) + E,f - mBo' 



The neutrino energy Ey is obtained from the measured value of Emiss (see Section [31] ). 
The agreement between data and simulation on Emiss has been verified using the side 
bands of the Ds£ sample. In order to have enough statistics to perform this test, cuts 
not correlated with the missing energy have been relaxed. In addition, to verify the 
resolution on the energy estimate, the studied sample has been enriched in light quark 
events by applying an anti-b-tagging cut (Section [2.41 ). A relative shift of A(MC — 
Data) = 500 MeV has been measured and the simulation has been corrected. Figure |^ 
shows the agreement between data and simulation after having applied this correction. 
The neutrino E^, resolution has been improved by correcting E^iss by a function of (Dgf) 
energy^ and which has been determined from simulated signal events: 

E, = Em^ss + F{E{^DJ)). 

The data-simulation agreement on p-^o has been verified on the selected signal events 
after subtraction of the combinatorial background (estimated from events selected in 
side-bands of Dg and signals) (Figure H). 



3.7 Proper time resolution and acceptance 

The predicted decay time distributions have been obtained by convoluting the theo- 
retical distributions with resolution functions evaluated from simulated events. Due to 
different resolutions on the decay length, different parametrizations of the proper time 
resolution have been used for three different classes in the Dg^ sample: KgK"*" decays, 
other non-leptonic decays and semileptonic Dg decays. Different parametrizations have 
been also used for the two Vertex Detector configurations installed in 91-93 and in 94-95. 

^here Dg means "observed decay products of Ds", including also the decays where the Ds is not fully reconstructed: 
specifically Da+ — > if>t^Vi and Ds"*" — > <j>h'^ X 
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The proper time resolution is obtained from the distribution of the difference between 
the generated (t) and the reconstructed (tj) time. The following distributions have been 
considered: 

• TZbiit — ti) is the resolution function for direct semileptonic B decays. 
TZbi is parametrized, for the Dgi sample, as the sum of three Gaussian distributions. 
The width of the third Gaussian is taken to be proportional to the width of the 
second Gaussian. 



- ti) = (1 - /2 - h )Git - U, ai) + hOt - U, a2) + hGit - U, as) 

m \ _2 i9 



with 0"! = J(y\i + cTpi^^ 
^2 = \jcrl2 + crl^t^ 

0"3 = 530-2 



In the (j)ih analysis a fourth Gaussian distribution has been added. 
The parameters related to the decay length and proper time resolutions, ai- and ap- 
respectively, and the relative fractions fi are listed in Table p. A typical parametriza- 
tion of the resolution, for the Ds£ sample, is shown in Figure |^ for the ^vr"*" decay 
mode obtained with the 94-95 Vertex-Detector configuration. 



IZhci is the resolution function applied to "cascade" events. 

Since the charm decay products have been only partially reconstructed in these 

events, the momentum of the B° candidate is underestimated giving a long positive 

tail in the proper time resolution function. 

The function, TZbciit — ti), is well described by a Gaussian distribution convoluted 

with an exponential distribution. The variation of the shape of this distribution 

with the generated proper time has been neglected. 



Dsi sample 


Dg decay mode 


(^Liips) 


o-pi 


(rL2ips) 


(^P2 


S3 


/2 


h 


K°K+ (92-93) 


0.16 


0.08 


1.04 


0.16 


- 


0.50 





KOK+ (94-95) 


0.16 


0.08 


0.98 


0.16 


- 


0.28 





other non-leptonic (92-93) 


0.11 


0.07 


0.39 


0.16 


5 


0.26 


0.07 


other non-leptonic (94-95) 


0.11 


0.07 


0.37 


0.16 


3 


0.16 


0.02 


(f)i+iy (92-93) 


0.14 


0.075 


0.31 


0.15 


6 


0.29 


0.09 


(f)i+iy (94-95) 


0.14 


0.075 


0.31 


0.15 


6 


0.21 


0.07 



(f)ih sample 


(^Liips) 


CTpi 


(7 mips) 


<7P2 


(^L^ips) 


0"P3 


(TLiips) 


crp4 


/i 


/2 


h 


0.13 


0.08 


0.28 


0.09 


0.32 


0.19 


1.06 


0.42 


0.31 


0.41 


0.10 



Table 8: Fitted values of the parameters of the resolution function TZbi{t — U) obtained, 
on simulated events, for the Dgi and (pih samples. 

Distortions on the reconstructed proper time can be due to a non-uniform reconstruc- 
tion efficiency as a function of the true proper time (acceptance). 
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Non-uniform efficiencies have been observed, on simulated events, in the Dg decay modes 
(f)TT, K*'^K and (piu because of the selection criteria on L/cr(L). 

This effect has been taken into account by inserting in the ffiting function, for those 
channels, an acceptance function {A{t)) parametrized on simulated events. 





Measurement of the B lifetime 



s 



The Bg meson lifetime has been studied using the signal sample (Section p.5|) and 
a background sample containing events selected in the sidebands of Dg (0) candidates. 
Sidebands events are "right" sign events lying in the Dg mass interval [1.91 — 1.93] U 
[2.01 — 2.15] GeV/c^ for the Dg hadronic decays and "right" sign events lying in the 
mass interval [0.990, 1.005] U [1.035, 1.060] GeV/c^ for the Dg semileptonic decays. 
In the Dg^ analysis "wrong" sign candidates have been also included in the background 
sample. This background sample is assumed to have the same proper time distribution 
as the combinatorial background in the signal sample. This assumption has been verified 
using the simulation. The probability density function used for events in the signal region 
is given by: 

-P(^i) = fbiPb/i'ti) + fbiPbli^i) + fbclPbciiU) + ffiPfi{ti) + fbkgPbkgiU). 

where ti and t are the measured and true proper times respectively. 
The different probability densities are expressed as convolutions of the physical probabil- 
ity densities with the appropriate resolution (TV) and acceptance (A) functions: 



for the signal: 



PuiU) = — expi-t/TBo)Ait) (g) TZuit - U) 



for the background coming from non strange B mesons: 



.B. 1 



Pm{U) = E fu — exp{-t/TBjA{t) ® TZuit - U 



q^s ^B. 



where q runs over the various B-hadrons species contributing to this background, 
for the "cascade" background: 

PbceiU) = J2 & — expi~t/TBjAit) ® TZbdit - U) 

for "fake lepton" candidates the function Pfiiti) has been parametrized using simu- 
lated events; 
for the combinatorial background two different parametrizations have been used: 

— Dg£ sample. 

PLgiU) = f-^exp{-t/r-) ® G{t - ti, a,) + f+^exp{-t/T+) ® G{t - U, a,-) + 
(l-r-/+)G(t-t„a,) 

Three distributions have been used for each of the three classes of decay time 
resolution aj (j = 1,3) (see Section |3?7|) . A negative exponential for poorly 
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measured events (with negative lifetime r~), an exponential distribution for the 
flying background (with lifetime r"*") and a central Gaussian for the non-flying 
one. The seven parameters (/~, /"'', r+, r^ and aj (j = 1,3)) have been fitted 
independently for the 92-93 and 94-95 data samples. The parameter aj (j = 1, 3) 
are taken to be different for the three classes of decay time resolution. 
— 4)ih sample. 

The combinatorial background shape has been described with a sum of four 
smeared exponentials {exp(ti, r) ® G{t — ti, cr)). 

Bg lifetime fit has been performed simultaneously on the signal and background samples. 
All parameters describing the shape of the background time distributions in the Dgi 
and (f)ih samples are left as free parameters. Results of the fit are shown in Figure [^ 
(Dsi sample) and in Figure ^ {(pih sample). Table ^ summarizes the different lifetimes 
measurements with their statistical errors. 



Decay mode 


Data set 


rB'iips) 


BJ 


Ds^07r 


(92-95) 


1 44+0.26 


Ds£ 


Ds ^ K''K+ 


(92-95) 


1 o 1+0.30 
-'-•<J-L-0.25 


BJ 


Ds ^ KgK+ 


(92-95) 


1.43l°il 


BJ 


Ds ^ K'''K'^+ 


(94-95) 


i.ooigj? 


DJ 


Ds -^ (pTT+rr^ 


(94-95) 


IMtni 


Ds£ 


Ds —>■ (f)7r~^7r^n~ 


(94-95) 


l.9Qtlil 


D,i 


Ds -^ (pi+u 


(92-95) 


-. 40+0.34 

-'-•'±^-0.27 


(pm 


(94-95) 


1.41 ±0.68 



Table 9: Bg lifetime determinations using the Dgi and (pih events samples. 



4.1 Systematic errors on the Bg lifetime 

Systematic uncertainties attached to the B^ lifetime determination are summarized in 



Table 10 



The main contributions to the systematic uncertainties come from: 
• Systematics from the evaluation of the B^ purity. 



Dsi sample: 

The different fractions for signal and background events have been calculated 
on an event by event basis. The expressions defining the effective purities are 
given in Section p.5.1| . The value of fhkg. has been varied according to the 



statistical uncertainties of the fitted combinatorial background fractions present 
in the different, Ds or K^K~, mass distributions. The value of fbd has been 
varied according to the errors given in Table ^ and in Table ^, which takes 
into account both the statistical error from the simulation and the errors on 
measured branching ratios. 
The evaluation of the systematics due to the procedure used to evaluate the 
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Systematics 


tbo variation in ps 


Jbkg. 

fbci 
fB 

Jbe 
Xd^ discrim. var. 

Pt discrim. var. 

rB+(1.65±0.04ps) 

rBo(1.56±0.04ps) 

t resolution 

t acceptance 

Simulated evts. statistics 


+0.0090 
-0.0130 
-0.0100 
+0.0110 
-0.0020 
+0.0020 

±0.008 
±0.004 

-0.0010 
+0.0010 
-0.0012 
+0.0013 

±0.008 
±0.010 
±0.020 


Total Syst. 


±0.03 



Table 10: Different contributions to the systematic uncertainty attached to the B^ lifetime 
measurement. 



Bg purity on a event by event basis has been evaluated in two steps. The 
distributions of the variable X^^ (Figure ^ for signal and background events 
have been re-weighted with as linear function in order to maximize the Data- 
simulation agreement: 



S{X^ 
B{Xj,^ 



B{Xn. 



bX 



Ds 



s 7 old 



The linear behaviour of the correction has been chosen because of the limited 
statistics in the data: it has been verified that a quadratic correction does not 
change the result significantly. 

The fit has been redone with this new probability distribution and the variation 
of the fitted lifetime value (±0.008 ps) has been taken as the systematic error. 
Because of the agreement between data and simulation (Figure |^e and ^f ) for 
the Pt distribution, the systematic error associated to this variable has been 
evaluated varying its distributions by the uncertainties of the parametrization 
obtained from simulated events. 
— (p£h sample: 

In this analysis the fractions of signal and background events have not been 
calculated on an event by event basis. The systematic uncertainty due to the 
variation of the fbce, fu ^^nd fhkg fractions have been obtained by varying these 
parameters by the errors reported in Table [^ and in Table ^. The systematic 
uncertainty attached to the fji fraction, affecting only the 0£ sample, has a 
negligible effect on the global result. 

Validation of the fitting procedure using simulated events. 

The fitting method has been verified on pure B^ simulated events: the measured 

value on this sample has been T-QQ{Dfl^y^^ = (1.605 ± 0.020)ps in agreement with 

the generated value [t^o = 1.6 ps). The statistical error of this verification has been 

included in the systematic uncertainties. 

A similar check has been performed on the (pih sample giving T^o{(j)ih.)^'~" = 
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(1.65 ± 0.04)ps. Since the statistical weight of the ipih channel is small compared to 
the full sample, the error on the fitting procedure is dominated by the statistics of 
Ds£ simulated events. 

Systematic from the proper time resolution. 

Uncertainties on the determination of the resolution on the proper time receive two 
contributions: one from errors on the decay distance evaluation and the other from 
errors on the measurement of the B^ momentum. The agreement between real and 
simulated events on the evaluation of the errors on the decay distance has been veri- 
fied by comparing the widths of the negative part of the flight distance distributions, 
for events which are depleted in B-hadrons. The difference between the two widths 
has been found to be of the order of 10%. 

The systematic on Bg momentum has been evaluated by comparing the momen- 
tum distribution on simulated events with the distribution, background subtracted. 



obtained from the data sample (see Section p.5.1| ). Effects from shift and width dif- 
ferences between the two distributions have been considered by changing the shape 
of the distribution of simulated events; it has been found that the main systematics 
comes from difference in width: the width on data has been estimated to be larger 
by a factor 1.07 ± 0.04. 

Taking into account these two effects the uncertainty on the time resolution has 
been, conservatively, evaluated by varying the parameters ai, and cxp. of the resolu- 
tion functions (see Table ||) by ±10%. 

Uncertainties on the acceptance determination have been also considered: the pa- 
rameters entering in the definition of the acceptance function have been varied ac- 
cording to the errors given by the fit on simulated events. 

The final result is: 

Tbo = lA2t^^\l{stat.) ± 0.03(syst.) ps. (8) 

5 Lifetime difference between Bg mass eigenstates 

The Bg (or B^) mesons are superpositions of the two mass eigenstates: 

|B°) = i=(|B^) + |B°)) ; |B°) = i=(|B°i) - |B°)). 

The probabihty density for N semileptonic Bg decays is proportional to: 

^ oc (Br(B° ^ ^X)rHe-^H* + Br(B° ^ iXWi^e-"^^') (9) 

at 

where MK^ ^ ^X)= r(BOj(L) - ^X)/rH(L). 

The semileptonic partial widths for B^ and Bl are assumed to be equal since only CP- 
eigenstates could generate a difference (semileptonic decays are not CP-eigenstates). 
It follows that the two exponentials are multiplied by the same factor and the probability 
density for the decay of a Bg or B^ at time t is given, after normalization, by: 

V{t) = -I^ (e-^H* + e-^^*) (10) 

J- H + J- L 

where Pl = Pb'' + APbo/2, Ph = Pro - APbo/2. 
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Two independent variables are then considered: r = I/FboQ and 6 = AT-qo/T^o. 
As the statistics in the sample is not sufficient to fit simultaneously r and 6, the method 
used to evaluate 6 consists in calculating the log-hkelihood for the time distribution 
measured with the Dsi and (pih samples and deriving the probability density function for S 
by constraining r to be equal to I/Fbo = tb,, = (1.56±0.04) ps fT^jQ (|rB,/rBd-l| < 0.01 
is predicted in [jl[). 

The log-likelihoods function described in Section ^ have been modified by replacing 
the physical function B^ (exp(— t/rBo)) by Equation ([T0|) and they have been added. 
The log-likelihood sum has been minimized in the (r, 6) plane and the difference with 
respect to its minimum (A£) has been calculated (Figure |TO|-a): 

AC = - log Clf^'^'\T, 6) + log ^^//■"^"((r)™", (5)™") . 

The probability density function for the variables r and 6 is then proportional to: 

V{t,6) oce"^^ 

The 6 probability distribution is obtained by convoluting V{t, 6) with the probability 
density function /(t-=t o )(''")' expressing the constraint r = t^o, and normalizing the 

d 

result: 

V{5) 



IT^i^,S)f{T=T^o)i^)drd5 



where 



/(r=r3o)(^) = l/(v^^rgo) exp(-(r " ^-B" ) VSc^r^o - 



The upper limit on AFBo/rBO, calculated from V{6), is: 

AFbo/Fbo < 0.45 at the 95% C.L. 

This limit takes into account both statistical uncertainties and the systematic coming 
from the uncertainty on the B^ lifetime. 

The systematic uncertainty originating from other sources has been evaluated by convo- 
luting V{t, 6) with the probability density functions of the corresponding parameters: 

J / [T, 0, Xgyg, ..., X gyg) J (j-=r-gp)\^ J J \^ sys) ■ ■ ■ J (^ sys)'^'^'^^ sys- ■ ■'^•^ sys 

V{6) = 

J / [T,0,Xgyg, ...^Xgygj j(r = T-gf))\j)j\XgygJ--.j[Xgygj(lTaXgyg...(lXgyg(lO 

d 

where x\yg are the n parameters considered in the systematic uncertainty and f{x\yg) are 

the corresponding probability densities. 

Since the method implies heavy numerical integrations over a n-dimensional grid only two 

systematics have been considered here: the purity in B^ meson of the selected sample 

and the acceptance. This approximation is justified since systematic uncertainties are 

expected to be small (as they are in the lifetime measurement) and dominated by these 

two parameters. 

The AFbo/Fbo probability distribution, obtained with the inclusion of the systematics, 

is shown in Figure 0-c, the most probable value for AFbo/Fbo is and the upper limit 

at 95% confidence level is: 

AFbo/Fbo < 0.46 at the 95% C.L. 



■* T does not coincide with the measured B^ lifetime if AFgo is different from zero 
* It has been assumed that AFb , = 0. 
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It should be noted that the world average of the B^ lifetime cannot be used as constraint 
in such analysis, since it depends on AFbo and on Fbo. Moreover, this dependence is also 
different for different decay channels. In the Dgi case the expression of the average B^ 
lifetime is given by: 

Tbo p^f = 'i ^ ^ 11 

rBo(i-(2ArBo/rBo)2) 
6 Study of Bg-B^ oscillations 

The study of B^-B^ oscillations requires the tagging of the sign of the b quark in the 
Bg meson at the decay and production times. 

The algorithm used for the b(b) tagging at production time has been tuned in order to 
have the best performances on the D^i sample, where all the charged particles from the 
Bg decays have been reconstructed. 



6.1 6(6) tagging at production time 

The signature of the initial production of a b(b) quark in the jet containing the B^ or B° 

candidate is determined using a combination of different variables which are sensitive to 

the initial quark state following the same technique as in Section p.l| . For each individual 

variable Xi, the probability density functions fb{Xi) (/^(Xj)) for b (b) quarks are obtained 

from the simulation and the ratio Ri = fj;{Xi)/fh{Xi) is computed. The combined tagging 

variable is defined as: 

1 p 

xtag = , , p , where R = Y[Ri. (12) 

i + Jx 

The variable Xtag varies between -1 and 1. High values of Xtag correspond to a high 
probability that a given hemisphere contains a b quark in the initial state. If some of the 
variables Xj are not defined in a given event, the corresponding ratios Ri are set to 1, 
corresponding to equal probabilities for the initial state to be b or b. 
An event is split into two hemispheres by the plane passing through the beam interaction 
point and perpendicular to the direction of the Bg candidate; then nine discriminant 
variables have been selected for this analysis. Five variables are defined in the hemisphere 
opposite to the Bg meson, in which reconstructed charged particles have been used: 

• the mean hemisphere charge which is defined as : 

In this expression n is total number of charged particles in the hemisphere, qi and 
Pi are, respectively, the charge and the momentum of particle i el is the unit vector 
along the thrust axis and k,=0.6; 

• the weighted sum of the charges of particles with tracks identified as kaon candi- 
dates: 



QK = Eqi{\p^■es\)''■, 

the sum of the charges of tracks having significant impact parameters with respect 
to the event primary vertex; 
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• the sum of the charges of the particles whose tracks are compatible with the event 
primary vertex; 

• the momentum transverse to the jet axis multiplied by the charge of the identified 
lepton candidate with the highest momentum. 

These variables have been combined to form the discriminant variable x" . 

Another set of three variables are evaluated in the hemisphere which contains the B^ 
meson candidate and only tracks not included in the B^ candidate decay products have 
been used in their determination Q. They are: 

• the mean hemisphere charge, computed using (|13D with el directed along the recon- 
structed momentum of the Bg candidate; 

• the rapidity with respect to the direction of the thrust axis multiplied by the charge 
of the identified kaon candidate with the highest momentum having a trajectory 
compatible with the primary vertex (this algorithm aims at reconstructing the frag- 
mentation kaon produced with the B^, this kaon has a sign opposite to the b quark 
contained in the meson); 

• the momentum of any reconstructed A'' candidate multiplied by the charge of the 
proton from its decay (same principle as in the previous item when a baryon instead 
of meson is produced). 

These variables have been combined to form the discriminant variable xf . In addition 
the distribution of the polar angle of the direction of the thrust axis, common to both 
hemispheres, is also used to benefit from the forward-backward asymmetry of the b quark 
production relative to the electron beam axis. 



6.2 Measurement of the tagging purity in events with an exclu- 
sively reconstructed D* 

The high statistics sample of exclusively reconstructed D*, accumulated in 1994-95, 
has been used to check the tagging procedure. The purity of the tagging at production 
time, etag, has been measured on those events using the analysis of the B[} — B|j mixing. 

The D*^ candidates have been selected by reconstructing the decay chain D*^ -^ D°7r+ 
followed by D° -^ K~7r"'" or D" -^ K~7r"'"7r°. The selection criteria rely mainly on the small 
mass difference between D*"*" and D° mesons [^. The measurement of the B|j — B^ mixing 
is performed by correlating a) the sign of the D* charge, which tags the B flavour at 
decay time (since D*~ in these events are mainly produced from B° and D*"^ from B 



with b) the global tagging variable, x° , evaluated in the hemisphere opposite to the D* 
and obtained by combining the five first quantities defined in the previous section. If 
the Bj meson, decaying into a D* , has oscillated, the D* charge and the value of the 
variable x°^g are expected to be of unlike sign. The mass difference, Am-go, between the 

two physical states of the B° — Bj system is obtained from the study of the D° decay 
distance distribution for unlike and like sign events. Details of the analysis can be found 
Q . The amplitude of the time dependent oscillation is sensitive to the probability of 



m 



correctly tagging events as unmixed and mixed B^J candidates. A fit has been performed, 
fixing the mass difference Am-QO to the world average [^, and leaving etag as a free 



parameter. The fit has been repeated for different minimum values of the global tagging 

^ In the Dsi? analysis all the Bj? decay products are identified and removed, for more inclusive analyses this is only 
partially possible 
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variable x° . Results are reported in Table |TT|, together with the predictions from the 



simulation. The fraction of events f events remaining after the cut on the tagging variable 
is also reported. 





Data 


Simulation 




^tag 


J events 


^tag 


J events 


Xia, >0.0 


0.68 ±0.02 


1.0 


0.69 


1.0 


<J>o.i 


0.69 ±0.02 


0.88 


0.71 


0.89 


^1a, >0.2 


0.71 ±0.02 


0.77 


0.74 


0.78 



Table 11: Values of etag obtained from the analysis of exclusively reconstructed D*^ for 
different cuts on the value of the tagging variable x° . Also reported is the fraction of 
events remaining after the cut. Expectations from the simulation are also given. 



The tagging efficiency, estimated using the D* sample, is consistent within its uncer- 
tainty with the expectation from the simulation. 

The selected sample of exclusively reconstructed D* still contains a significant fraction 
of events originating from charm and light quarks. In order to study the distribution of 
the tagging variable x^^^, the b-tag probabihty for all tracks in the event has been required 
to be smaller than 10~^ |T^. The fraction of non-b events in the remaining sample is 
estimated to be 5%. The distribution of the product between the D*^ charge and the 
value of the tagging variable x° is shown in Figure |ll]-a, together with the expectation 
from the simulation. 

Another check has been performed by selecting events with an exclusively recon- 
structed D* accompanied by a lepton of opposite charge. This sample is highly enriched 
in Bj, but has a limited statistics. However, it allows the study of the tagging variable 
xl defined, in the same hemisphere as the D* -lepton candidate, by combining the other 
three variables mentioned in the previous section. The variable which quantifies the pres- 
ence of an identified kaon of highest momentum compatible with the primary vertex has 
been removed from the definition of xl . The distribution of the product between the 
D* charge and the value of the tagging variable, x^ , is shown in Figure |TT|-b together 
with the expectation from the simulation. 

The selected Ds£ sample do not have enough statistics to perform a quantitative check. 
The Xtag distributions expected from the simulation and measured data, using the Ds£ 
sample, are found to be compatible within statistics (Figure 0). 



6.3 Tagging procedure 

An event is classified as a mixed or an unmixed candidate according to the relative 
signs of the Dg electric charge, Qd, and of the Xtag variable. Mixed candidates have 
Xtag X Qd < 0, and unmixed ones x^g x Qj) > 0. 

The probability, e?,, of tagging the b or the b quark correctly from the measurement of 
Xtag has been evaluated using a dedicated simulated event sample and has been found to 
be, in the Ds£ sample, 74.5 ± 0.5% in 94-95 data and 71.5 ± 1.2% in 92-93 data. 
In the (f)lh sample the tracks from the B decay have not been all reconstructed. The 
tagging purity is lower with respect to the one estimated in the Dg^ sample due to some 
possible misidentification between primary and secondary tracks present in the same 
hemisphere as the (j) meson. The value found in simulated events is (e^ = 0.69 ± 0.01). 
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To improve the tagging purity further, the shape of the Xtag distribution can be in- 
cluded in the analysis. 
Four purities enter in the analysis: 

- Em'- tagging purity for the direct b ^ i decays; 

- ebci'- tagging purity for 6 ^ c — *> £ "cascade" decays; 

" e,!^j^g "^^^ probabihty of classifying background candidates as mixed or as unmixed 
(computed on sidebands events); 

- ^^^yunrmx) Y)xo\:,ahiYitY of classifying fake lepton candidates as mixed or as unmixed. 

using Xtag as a discriminant variable each of these purities is replaced by the function 

eX{xtag)^ where X is the Xtag probability density function. 

The global probability density function has been divided by the sum eXl^{xtag) + (1 — 

e)XJ^^{xtag) {r = right tag and w = wrong tag) in order to keep, for the signal part, the 

relation e"" = 1 — e'^' . 

The functions entering in the final likelihood are then re-defined as: 

-vr _ ^bl^Mv^tag) y'w _ i y 



^M^ui^tag) + (1 ~ ^bd^ui^t, 



b£ — ^ y^be 



agj 



vr _ Cfcrf^fcrf vw _ (1 £brf)^fce£ 

^''-^ e,eXUxtag) + (1 - eb,)X-(xi,3) ^^-^ euX^Xtag) + (1 - eM)X^Axtag) 

,mix vmix ^unmix vunvaix 
Ymix __ ^bkg bkg ^1^1^"^^^ = ^fcfcg bkg 

^o'^a euXl,{xtag) + (1 - eu)X^,{xt^,) ^^^^ euX^Xtag) + (1 - e,,)X-(xt,,) 

^mix ymix ^unmix yunmix 
j^mix ^ ^fi ^fi j^unmix ^ ^/^ ^fi 

■' ^be^bei^tag) + (1 "~ ^bd^ui^tag) ^ ^bi^bii^tag) + (1 ~ ^b£)Xbi{Xtag) 

The effective tagging purities obtained, in the Dgi sample, with this method correspond 
to 78 ± 0.5% for 94-95 data and to 74 ± 1.2% for 92-93 data. 

6.4 Fitting procedure 

From the expected proper time distributions and the tagging probabilities, the prob- 
ability functions for mixed and unmixed events candidates have been computed Q 

^™^(^.) = fSPSo^t,) + ff,Pr{td + AcATl^.) + ff£PTr{t^) + fbkgP^g^U). (14) 

where tj is the reconstructed proper time. The analytical probability densities are as 
follows, with t being the true proper time: 

• Bg mixing probability. 

PBfiU) ={ Xl,V^f{t)+X^,V^r'''{t) M(t)®7^5,(i-i^) (15) 

• "cascade" background mixing probability. 



(/fJ/2)( xi^v^p'^it) + x]^^,r^fit) )+ 

(/|/2)( X^^,V~%t) + Xl^V^fit) )+ (16) 

fb^^bJ^B+ exp(-t/rB+) 

ft^XljT^^ exp(-t/rAj }A{t) ® Ubciit - U) 



T 



In the following, only the probability function for mixed events is written explicitly; the corresponding probability for 
unmixed events can be obtained by changing r — » ui. 
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Note that the two terms contributing to the Bg are due to the fact that, in the decay 
Bg -^ D+D~X, the lepton can originate from either Dg mesons. The Bg contribution 
can then be simphfied in the expression and becomes, mixing independent: 

/Bo/i"B2exp(-t/rBo) 

• non-strange B-hadrons mixing probabihty. 

/rX,-,/rB+exp(-t/rB+) (17) 

fttXfjT^, exp(-t/rAj }A{t) ® nu{t - U) 

• mixing probabihty for candidates from hght quark events or fake leptons: 

Vjr{ti) = XfrVf,{U) (18) 

• combinatorial background mixing probabihty: 

^bkgi'ti) = ^bkg'^bkgiti) (19) 

The parameters entering in the proper time distribution for this background have 
been determined in the hfetime fit. 

The osciUation analysis has been performed in the framework of the amplitude method 



25| which consists in measuring, for each value of the frequency Attibo , an amplitude A 
and its error <j{A). The parameter A is introduced in the time evolution of pure B^ or 
B^ states so that the value A= \ corresponds to a genuine signal for oscillation: 

P(B° -^ (B°, bJ)) = —e-^s X (1 ± A cos(AmBot)) 

The 95% C.L. excluded region for AmBO is obtained by evaluating the probability that, 
in at most 5% of the cases, a real signal having an amplitude equal to unity would give an 
observed amplitude smaller than the one measured. This corresponds to the condition: 

yl(AmBo) + 1.645 or(A(AmBo)) < 1. 



In the amplitude approach it is possible to define the exclusion probability, that is the 
probability that a certain Am-QO value lies in an excluded region if the generated AmBO 
was very large (ArraBO — * cxd). The sensitivity is the value of Am-BO corresponding to 50 
% of exclusion probabihty. 

Using the amplitude approach (Figure |13D, and considering only statistical uncertain- 



ties, a limit has been obtained: 

AmBO > 7.4 ps-i at 95% C.L. (20) 

with a corresponding sensitivity at Att^bo = 8.3 ps^^. At Att^bo = 10 ps~^, the error on 
the amplitude is 0.85. 

Several checks have been done to verify the reliability of the amplitude fit: the proper 
time distributions for mixed and unmixed events have been verified to be well reproduced 



by the fit (Figure |Tj-a,-b) and the ratio between mixed events and the total number of 



events in bins of the proper time has been compared with the expected distribution for 
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ArriBO = 5 ps'^ and Am^o = 10 ps~^. These values have been chosen to illustrate the 
behaviour of the expected oscillation curve for A = {Am^o = 5 ps~^) and A = 1 
(A?7iB0 = 10 ps~^) (Figure |l^e). It could be seen that the oscillation curve at AmBO = 
10 ps~^ (where A is close to 1) fits the data better than the corresponding curve at 
Am-BO = 5 ps~^ (where A is compatible with 0), as expected from the definition of A. 

6.5 Study of systematic uncertainties 

Systematic uncertainties have been evaluated by varying the parameters which have 
been kept constant in the fit, according to their measured or expected errors, using the 
formula 



a[Aii.)]=AAiu) + il-A)-^''^'^'^ 



a[Aiu)y 



AA{u) and Aa{A) indicate the variations of the amplitude, in the central value and in 

the error, due to the considered systematics. 

Three main sources of systematic uncertainties have been identified: 

• Systematics from the tagging purity. 
— Dgf sample. 



The studies done in Section [63 show that, using the tagging variables in 
the opposite hemisphere and requiring \x°^g\ > 0., the difference between 
the values of the tagging purity measured in real and simulated events is 
etag{DATA) - etag{MC) = -0.01 ± 0.02. It has been verified that the real and 
the simulated distributions for the tagging purities agree in both hemispheres. 
The systematics coming from the control of the tagging purity has been evalu- 
ated by varying the probability distributions of the discriminant variable for b 
and b quarks in a way to induce an absolute variation on the effective value of 
the tagging purity of ±3.0%. 

— (j)£h sample. The agreement between data and simulation has not been checked 
for this sample; a conservative absolute variation of 5% in the tagging purity 
has been assumed. 

• Systematics from the 3° purity. 

The same procedure already applied for the lifetime measurement has been used. 

• Systematics from the resolution on the B decay proper time. 

The same procedure already applied for the lifetime measurement has been used. In 
addition, the systematic error due to the variation of the proper time distribution of 
the combinatorial background, has been considered: the parameters used to define 
the background shape, in the lifetime fit, have been varied according to their fitted 
errors. 

The inclusion of systematic uncertainties lowers the sensitivity to 8.1 ps~^ without 
affecting the 95% C.L hmit. In Table |T^the amplitude values are reported, together with 
their statistical and systematical errors, for five different values of Att^bo. 

The exclusion probability of AmBO = 7.4 ps~^ is 54% while the probability of obtaining 
a limit on Am-QO higher than the actual one is 38% (Figure p!5|-c). 



Figure |15|-a and Figure |15|-b represent, respectively, the error on the amplitude and the 



exclusion probability as a function of AmBO . 
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Am^o{ps ^) 


A 


aA{stat) 


cta (total syst. (but t resolution)) 


cr^(t resolution syst.) 


2.5 


-0.638 


0.304 


0.112 


0.033 


5 


0.037 


0.400 


0.118 


0.060 


7.5 


0.182 


0.561 


0.069 


0.098 


10 


1.343 


0.846 


0.160 


0.180 


12.5 


0.867 


1.241 


0.285 


0.389 



Table 12: Amplitude values with statistical and systematic errors for three different values 
of Athbo 

7 Conclusion 

A sample of 436 T)f£^ candidate events has been selected from about 3.6 million 
hadronic Z^ decays accumulated by DELPHI between 1992 and 1995, using seven different 
Dg decay modes. The number of events coming from B^ semileptonic decays has been 
estimated to be 230 ± 18 in this sample. In addition, a sample of 441 (pih, containing 
41 ± 12 Bg semileptonic decays, has been also used. Events contained in the Dsi sample, 
with a reconstructed and have been removed from this last sample. 
Using these samples, three analyses have been performed. The B^ lifetime has been 
measured and a limit on the fractional width difference between the two physical B^ 
states has been set: 

r(B°) = (1.42+°:i^(stat.) ± 0.03(syst.)) ps 

ArBo/EBO < 0.46 at the 95% C.L. 

This last result has been obtained under the hypothesis that t^o = tb^. 

The study of B^ — B° oscillations sets a limit at 95% C.L. on the mass difference between 

the physical Bg states: 



AmBO > 7.4 ps-i at 95% C.L. 



(21) 



with a corresponding sensitivity equal to 8.1 ps~^. 

Previous DELPHI results obtained with Dg^ and <pi samples (|]26|,[^) are superseded by 

the analyses presented in this paper. 
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Figure 1: The plot on the top shows the distribution of the X-^,^ discriminant variable for 

the (pTT channel in 94-95 data. The points with error bars represent the data, the white 

histogram shows the contribution from the simulated signal and the shaded histogram 

shows the contribution coming from simulated background events. 

It could be seen that the Xd^ is able to discriminate the signal (Dg) from the combinatorial 

background. 

The four figures on the bottom show the effect, on the 07r signal in the 94-95 data, of 

a cut on the discriminant variable (white histograms represent "right-sign" events while 

shaded histogram show "wrong-sign" events). 
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Figure 2: Invariant mass distributions for Dg candidates in six non-leptonic decay modes 



!)7r^ 



r°K- 



F*0 



K^K"*", K K*"*", (pTT^n^Ti^ and (f)7[~^TT^). The last three decay modes have 
been reconstructed using only the 94-95 statistics. The corresponding distribution for 
wrong-sign combinations are given by the shaded histograms 
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Figure 3: a) Invariant mass distributions for Dg candidates in non-leptonic decay modes 
((pn'^ , K K+, K°K+, K K*+, (pTT^ii^Ti^ and (pn^TT^). b) K^K~ invariant mass distribu- 
tion for Ds candidates selected in the two semileptonic decay modes. The corresponding 
distribution for wrong-sign combinations are given by the shaded histograms. The curves 
show the result of fits described in the text. 
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Figure 4: K^K^ invariant mass distribution of the (pih candidates. The curves show 
the result of the fit described in the text (the signal and the combinatorial background 
components are represented by the white and shaded histograms respectively). 
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Figure 5: The plots illustrate the agreement between data and simulation for the vari- 
ables used in the estimate of the B^ purity on an event by event basis. 
a),b),c) and d) show the X-q^ distributions for the channels (pii, K°*K, KgK and (f)i~^ 
respectively. White histograms and shaded histograms represent the signal and the back- 
ground respectively. 

e) and f) show the px distributions for the samples selected with pr > I GeV/c ((J)tt'^ and 
'v) and pt > 1-2 GeV/c (all the others) respectively. 



For the px distribution the Ds-D and the combinatorial background are considered sepa- 
rately. 
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Figure 6: The plot on the top shows the comparison between data and the simulation for 
the missing energy distribution after correction (see Section p.(\ ). The data (dots with 
error bars) and the simulation (dashed histogram) are enriched in light quark events. The 
dotted histogram shows the missing energy distribution in simulated b semileptonic decays. 
The plot on the bottom shows the comparison between the Bg momentum distribution for 
simulated events and the one estimated from data in the signal region by subtracting the 
Bg momentum distribution of events in the Dg side bands from that of the events in the 
signal region. 
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Figure 7: T/ie six plots on the top show the proper time resolution, defined as the dif- 
ference between the generated time (tsim) o-nd the reconstructed time (tree), in bins of 
generated time on a sample of (J)tt events simulated with the 94-95 Vertex Detector con- 
figuration. The plot on the bottom shows the distribution of the reconstructed time with 
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Figure 8: Dsi sample, a) Likelihood fit for events in the signal mass region. The points 
show the data and the curves correspond to the different contributions to the selected 
events, b) The same as a) but for "wrong-sign" events and for events situated in the side 
band region. 
The labels of the signal components refers to Section \3.5.^ 
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Figure 9: (j)ih sample. 

Likelihood fit for events in the signal mass region. The points show the data and the 

curves correspond to the different contributions to the selected events. 



The labels of the signal components refers to Section \3.5.i 
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Figure 10: a) 68%, 95% and 99% C.L. contours of the negative log-likelihood in the 
plane r(= l/rBj-AFBo/rBO evaluated on the Ds£ sample, b) Same as a) but with the 
constraint t = Tgo. c) Probability density distribution for AT^o/T-qo; the three shaded 
regions show the limits at 68 %, 95% and 99% C.L. respectively. 
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Figure 11: Check of the flavour tagging on the D* sample. 

a) Distribution of the product between the global tagging variable x° and the charge of 
the D* in the hemisphere opposite to the D* candidate. 

b) Distribution of the product between the global tagging variable xf and the charge of 
the D* in the same hemisphere as the D* -lepton candidate. 

The full dots with the error bars represent the data. The histogram is obtained in the 

simulation. 

The non perfect separation is due to the mistag fraction of Xtag but also to the B[J-B[j 

mixing. 
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Figure 12: The plot shows the distribution of the Xtag discriminant variable in the Ds£ 
sample. The points with the error bars represent the data, the white histogram shows the 
contribution from combinatorial background, the lighter histogram the contribution from 
DgD events and the darker histograms the contribution from the Bg signal in which Bg 
mesons produced from borb quarks have been distinguished. 
The degree of separation between the b and b histograms quantify the tagging purity of 

•^tag- 
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Figure 13: Variation of the oscillation amplitude A as a function of ArUs- The lower 
continuous line corresponds to A + 1.645 a^ where a^ includes statistical uncertainties 
only, while the shaded area includes the contribution from systematics. The dashed-dotted 
line corresponds to the sensitivity curve. The lines at A=0 and A=l are also given. 
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Figure 14: Proper time distribution of mixed and unmixed events in Dgf sample (a-h) and 
in the (pi sample (c-d); the full dots with error bars represent the data, the curves are the 
corresponding distributions for Am^o = 10 ps~^ . 

c) Ratio between the mixed events and the total number of events in bins of proper time in 
the Ds£ sample. The full (dashed) line represents the prediction for an oscillation (A = 1) 
with Am-BO = 5 ps^^{10 ps~^). 
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Figure 15: a) Variation of the error on the amplitude as a function of Arris ■ b) Exclusion 
probability vs. Anig. c) Lower limit probability density function vs. Anig. 



